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不 确定 混沌 系统 自 适应 改进 投影 同步 与 参数 估计 * 


FRE, R fi 
(西北 工业 大 学 理学 院 应 用 数学 系 ， 西 安 710072) 


摘 ”要 : 研究 一 类 新 混沌 系统 与 Genesio-Tesi 混沌 系统 的 改进 投影 同步 ， 并 对 新 混沌 系统 未 知 参 数 进行 估 
计 。 依 据 主动 控制 原理 ， 设 计 了 非 线 性 控制 器 和 参数 白 适 应 更 新 律 ， 用 Lyapunov 稳定 性 理论 证 
明了 误差 系统 原点 的 渐进 稳定 性 。 数 值 仿真 结果 验证 了 设计 的 非 线性 控制 器 和 参数 自 适应 更 新 律 
的 有 效 性 及 其 在 有 界 噪声 干扰 下 的 鲁 棒 性 。 
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混沌 是 非 线性 系统 的 特有 现象 ， 表 现 为 系统 长 期 演化 结果 敏感 依赖 于 系统 初始 条 件 。 基 于 
混沌 同步 在 保密 通信 等 领域 的 广泛 潜在 应 用 价值 ， 以 及 人 们 对 其 奇异 特性 的 探索 兴趣 ， 自 上 
世纪 九 十 年 代 Pecora 等 目的 开创 性 工作 以 来 ， 众 多 学 者 对 混沌 同步 展开 了 广泛 的 研究 ， 提 出 
了 混沌 同步 控制 的 多 种 方法 和 类 型 中。 常见 的 混沌 同步 主要 有 完全 同步 由、 相同 步 同 、 延 迟 同 
p, RAFE, JAHE, BEART 等 。 同 一 系统 的 混沌 同步 较 易 实现 ， 而 不 同系 统 
的 混沌 同步 则 较 难 。 关 于 随机 力 对 混沌 同步 的 影响 ， 部 分 学 者 大 多 使 用 高 斯 白 噪 声 进 行 过 研 
究 ， 然 而 实际 系统 所 承受 的 环境 噪声 大 多 是 有 界 的 ol 。 本 文 将 基于 主动 控制 原理 ， 设 计 非 线 
性 控制 器 和 参数 更 新 律 ， 酉 究 一 类 新 混沌 系统 A Genesio-Tesi 混沌 系统 四 的 改进 投影 同 
步 ， 估 计 新 混沌 系统 的 未 知 参数 ， 并 研究 有 界 噪 声 对 混沌 同步 和 参数 估计 结果 所 产生 的 影响 。 


2 模型 及 理论 分 析 
本 文 研究 的 新 混沌 系统 模型 方程 为 (14 
tı —Z273 十 QZ1， 
$2 = Z173 + br2, (1) 
t3 = iziz 十 cZ3. 
ZE X = (zxi,22, 23) 表示 系统 的 状态 ，a, b, c 是 系统 参数 。 在 a = 5, b= 一 10, c = 一 3.8 时 系 
统 (1) 有 图 1 所 示 混 沌 吸引 子 。 
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Genesio-Tesi 泥 沌 系统 方程 可 号 为 站 


Yı = Y2 + U1, 
fo = ys + U2, (2) 
ja = —a1yı — biy2 — C1Yy3 + Y? + ua. 

iX LY = (yi, y2,y3) FUU = (u1, u2, ua) 439] e o ZR BEGSCRDS] RAMAR a, bi, 


co 是 系统 参数 ， 满 足 条 件 Dhicl «a. Æa = 6, b = 2, c = L25IU =0( 即 对 Genesio- 
Tesi 系统 不 施加 控制 ) 时 ， 系 统 (2) 有 图 2 所 汞 混沌 吸引 子 。 


-20 -2 
x() 40 9  — yq) y 1 -* yi 


图 1: a—5,b—-10, c= —3.8 if, H2: a =6, bı 29, c1 = 1.2, 
系统 (1) 的 奇怪 吸引 子 系统 (2) 的 奇怪 吸引 子 


设 系 统 (1) 为 驱动 系统 ， 共 三 个 参数 a, b, ce 未知， 系统 (2) 为 响应 系统 。 下 面 研 究 系 
AOE j 系 统 (2) 的 改进 投影 同步 和 对 系统 (1) 三 个 未 知 参数 a, b, c 的 自 适应 估计 问题 。 
定义 1 Roa, oz Q3 为 三 个 不 同 的 尺度 因子 ， 当 下 面 三 个 等 式 同时 成 立时 ， 称 系统 (1) 与 
系统 (2) 实现 改进 投影 同步 [9 。 


lim lai) oun) =0, Jim, | a2() ~ oa) |= 0 jim |l aC) — ast [= 0. (3) 
G1 Ha = as = as =1 时 ， 称 系统 (1) 与 系统 (2) 实现 完全 同步 。 
注 2 当 oa = as = as = -1 时 ， 称 系统 (1) 与 系统 (2) 实现 反 同 步 。 
由 定义 1 知 ， 在 改进 投影 同步 意义 下 ， 系 统 (1) 与 系统 (2) 的 误差 动力 系统 为 
éi(t) = 一 z273 十 QZ1 — Q1Y2 — 013, 
é2(t) = Z173 十 bra — Q2Y3 — Q212， (4) 
é3(t) = iriza + cz3 + &3a1yı + a3 有 ge + Q3c1Y3 — a3y1 — 03t3- 
这 里 
ev) = ai (t) — ayı (t), e2(t) = zx2(t) ~ aayz(t), eslt) = zs(t) — osus(t), (5) 
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分 别 表示 系统 (1) 与 系统 (2) 各 状态 变量 间 的 同步 误差 。 从 而 系统 (1) 与 系统 (2) 的 改进 投影 同 
步 问 题 可 转化 为 误差 系统 (4) 的 原点 渐进 稳定 性 问题 。 这 里 u, ua, ua 为 待定 非 线 性 控制 函数 。 
为 在 实现 系统 (1) 与 系统 (2) 改进 投影 同步 的 同时 ， 通 过 系统 (2) 的 三 个 参数 al, bi, cl 的 自 适应 
变化 分 别 估计 出 系统 (1) 的 三 个 未 知 参数 o, bc， 依 据 主动 控制 原理 ， 非 线性 控制 器 可 取 为 
u(t) = d-[key — z2£3 一 aag2 十 a171]， 
uj(t) = a [kez + Z173 — Q2Y3 + 51232], (6) 


ua(t) — 去 [kses 十 iriz 十 asaly1 + 03b1Y2 + &3c1y3 ~ 03Yy? + eiza]. 


ZE k; > 0, i= 1,2,3。 参 数 自 适应 更 新 律 取 如 下 形式 


ài = Z1el， 
b, = Z262， (7) 
C1 = TZ3€3. 


3|38 1(Barbalat 引 理 )03-19 it f(t) 是 连续 光滑 的 实 函 数 ， 对 于 范 数 L2。 和 Lo。， 若 f(t) € 
L2 N Lo; 并 且 


则 有 lim. f(t) — 0. 

定理 1 对 于 系统 (1) 与 系统 (2) 的 任意 初始 条 件 和 任意 尺度 因子 01, az, as， 在 非 线性 控制 
器 (6) 和 参数 更 新 律 (7) 作用 下 ， 系 统 (0) 与 系统 (2) 能 够 实现 改进 投影 同步 。 

证 明 HW Lyapunov 函数 为 


V= 5$ e + e +e? +e + e2). (8) 
这 里 eu = al — a, ep = bi — b, ec = C1 一 Cc 分 别 表 示 用 a1, bi, ci 估计 未 知 参 数 a, b,c 时 的 误差 ， 
Mi éa = ài, éb = à, èc = 所。 利用 式 (6) 和 式 (7) 可 求 得 Lyapunov 函数 沿 误差 系统 (4) 的 导 
数 为 


dV 


rp €1€1 + e262 + £3€3 + eaea + ebeb + ecec 


(一 Z273 + azx1 一 o1y2 一 &1U1)e1 + (z123 十 br2 一 QA2Y3 一 02u2)e» 


ll 


十 (52s 十 cz3 + 030131 十 osbiya + ascia — osyi — asus) e3 
+å: (aı — a) + bı (bi — b) + &(ei — c) 

= (—khiei + azı — a121)e1 + (—Kae3 + 63 — 6122)e2 + (—k3e3 + x3 — c123)ea 
十 Ziel(al — a) + z2€2(b1 — b) + zaes(c1 — c) 

= —kye2 — kae2 — kze? = —e" Ke. 


这 里 e = (ei, e2, e3)T, K = diag(ki, ko, ka3)。 因 VV 是 半 负 定 的 ， 我 们 还 不 能 立即 得 知 误差 系 
统 (4) 的 原点 是 渐进 稳定 的 。 由 V < 0 Xll e, €2, €3, Ca, €b, €e € Loos 由 误差 系统 (4) 知 ĉi, €2, €3 
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€ Lo. WV --eTKex0, KEX, 从 而 
t i t , 
f Amin (K)||el?dt < f eTKedt < i -V dt = V (0) — V(t) < V (0). 
0 0 0 


XX HAS (X) XE XB PE K WEIER, Wei, e2,e3 € L2 B5114, Æt 一 oo 时 ， 
A e(t), ex(t), et) 一 0。 从 而 系统 (1) 与 系统 (2) 实现 改进 投影 同步 。 同 步 速 度 可 通过 1, ko， 
ka 进行 调节 。 
3 ”数值 仿真 

为 验证 上 文 设计 的 非 线 性 控制 器 (6) 和 参数 白 适 应 更 新 律 (7) 的 正确 性 ， 用 四 阶 Runge- 
Kutta 方 法 同时 求解 系统 (1) 与 系统 (2)， 计 算 步 长 取 0.001。 系 统 (1) 的 未 知 参 数 设 定 为 a = 
5, b = —10, c = -3.8， 初 始 条 件 zi(0) = —1, z2(0) = —1, z3(0) = -1。 系 统 (2) 的 参数 
取 al = 6, b; = 2.9, c = 1.2， 初 始 条 件 y(0) = 1, y2(0) = -3，%a(0) = 2。 尺 度 因 子 
Wa = —2, az = 0.6, a3 = —3. RA (1) 与 系统 (2) FG RAPRA (5) 计算 。 同 步 速 
ETRA ki = 15, ko = 10, ka — 2。 先 让 系统 (1) 与 系统 (2) 自由 演化 10 秒 以 使 其 达到 充 
分 混沌 状态 ， 在 上 = 10 秒 时 打开 控制 器 开关 ， 两 系统 状态 变量 的 同步 误差 变化 如 图 3 所 示 ， 
对 参数 a, b, ce 的 估计 结果 如 图 4 所 示 。 由 图 3 知 ， 在 上 = 14 秒 左右 时 ，e1(t) 5 e(t) 就 已 经 收敛 
到 零 ，es(t) 则 收 化 到 零点 很 小 的 邻 域 内 。 在 t = 24 秒 左右 时 ，es 人 ( 也 准确 收敛 到 零 ， 从 而 系 
统 (1) 与 系统 (2) 严 格 实现 改进 投影 同步 。 由 图 4 知 ， 在 上 = 16 秒 左右 时 ， 系 统 (1) 的 三 个 未 知 
2 a, b, c 就 被 ai (办 (办 ，cl( 轴 通过 白 适 应 变化 准确 估计 出 来 。 


et 


a(t), b,(t), c,(t) 


5 10 15 20 25 30 35 40 
Us 
E T 
7905 5 10 15 20 25 36 35 40 
t/s 
图 3: 系统 (1) 与 系统 (2) 的 改进 投影 同步 误差 图 4: 系统 (1) 三 个 未 知 参数 的 估计 值 


4 “有 界 噪声 影响 


由 于 系统 在 实际 环境 中 难免 要 受到 噪声 的 干扰 ， 而 这 些 噪声 往往 是 有 界 的 ， 因 而 研究 有 界 
噪声 对 系统 (1) 与 系统 (2) 改进 投影 同步 和 参数 估计 的 影响 更 符合 实际 情况 。 有 界 噪声 模型 可 取 
HE) = Dcos(Qt-- ou 人 十 FT)， 这 里 万 为 有 界 噪声 强度 ，Q 表示 激励 的 中 心 频率 ，uw 人 t) 是 标 
准 Wiener 过 程 ，o 表示 频率 随机 扰动 强度 ，T 是 服从 fo, 22] 上 均匀 分 布 的 随机 相位 15-171。 系 
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统 (1) 与 系统 (2) 在 有 界 噪声 &(t) 激励 下 的 模型 方程 可 分 别 表 示 为 
4, = —r273 + axı + Ẹ(t), 
Xo = T1723 + br2 + E(t), (9) 


£3 = iriz2 T cz3 + E(t); 


jı = yo + ui + K(t), 
yz = ys + uz + E(t), d) 
ys = —aiyi — biya — cys +y? + us + E(t). 


由 于 在 噪声 干扰 下 ， 系 统 (9) 与 系统 (10) 的 输出 可 看 作 是 随机 过 程 ， 央 此 ， 
两 系统 的 状态 变量 误 状 el (t), e2(t), ea(t) 严格 等 于 零 往往 很 难 实现 。 从 实用 角度 考虑 ， 


m. lex (t)] <7, k = 1,2,3, 


时 即 可 认为 系统 (9) 与 系统 (10) 实现 了 同步 (这 里 7 > 0AE, XO IRIZPAR AU). F 
面 研究 有 界 噪声 对 系统 (9) 与 系统 (10) 改进 投影 同步 ， 以 及 对 未 知 参 数 白 适应 估计 值 的 影响 。 

取 有 界 噪声 强度 D = 1.2, Q = o = 1 (其它 参 数 及 初始 条 件 同 前 )。 先 使 系统 (9) 与 系统 (10) Á 
由 演化 10 秒 ， 在 上 = 10 秒 时 打开 控制 器 开关 ， 两 系统 的 状态 同步 误差 变化 如 图 5 所 水 。 对 参 
数 a, b, ce 的 估计 结果 如 图 6 所 示 。 由 图 5 知 ， 在 上 = 14 秒 左右 时 ， 系 统 (9) 与 系统 (10) 三 个 状 
态 变量 的 状态 同步 误差 均 已 收敛 到 零点 附近 很 小 的 邻 域内 ， 同 步 误 差 曲 线 略 显 粗糙 ， 表 明 系 
统 (10) 的 状态 变量 可 以 按照 给 定 的 尺度 因子 ， 以 较 高 的 精度 跟踪 上 系统 (9) 相应 状态 变量 的 变 
化 。 跟 踪 误 差 是 一 个 非常 小 的 量 。 由 图 6 可 以 看 出 ， 在 上 = 15 秒 左右 时 ， 未 知 参 数 a, b, cC 
值 即 被 ul (£), 5 人 ,ci 人 通过 自 适 应 变化 以 很 高 的 精度 估计 出 米 ， 三 条 估计 值 曲线 略 显 粗糙 ， 

估计 误差 是 一 个 很 小 的 量 。 为 定量 刻画 有 界 噪 声 对 系统 (9) 与 系统 (10) 同步 误差 和 参数 估计 相 
对 误差 的 影响 ， 记 


absek = max |zk(t — oxyk(t), k= 1,2,3, 


分 别 表示 系统 (9) 与 系统 (10) 各 状态 变量 问 稳 态 时 的 最 大 同步 误差 。 


a = pag (|E x 100%), ea = mag {|= o. 


E3 = max {|2 =t] x 100%}, 
t>35 c t 


PARIRI RMBH a, b, cli VE E REDE TR 28 CE TS TE IRA HR «SETA FD EE UR S 
强度 D， 在 t = 10 秒 时 打开 控制 器 开关 ， 系 统 (9) 与 (10) 三 状态 变 景 问 稳 态 最 大 同步 误差 
值 和 对 未 知 参数 估计 值 的 相对 误差 的 最 大 值 如 表 1 所 示 。 由 表 1 可 见 ， 随 着 有 办 噪声 强 
度 刀 的 增 大 ， 最 大 同步 误差 和 参数 估计 相对 误差 的 最 大 值 呈 逐渐 增 大 趋势 。 在 有 界 噪声 强 
度 D < 6 时 ， 状 态 变 量 最 大 跟踪 误差 不 超过 0.9463， 这 相 比 于 状态 变量 的 变化 范围 是 非常 小 的 
一 个 量 ， 而 参数 估计 相对 误差 的 最 大 值 不 超过 3.7263%。 注 意 到 表 1 中 所 列 数据 是 系统 (9) 与 
系统 (10) 在 达到 稳 态 同步 时 同步 误差 和 参数 估计 相对 误差 的 最 坏 结果 ， 因 此 ， 若 用 稳 态 同步 
时 a1(t), b1(t), ex (t) 的 均值 或 中 位 数值 作为 对 未 知 参 数 a, b, c 的 估计 ， 其 相对 误差 将 会 显著 低 
于 表 1 中 所 列 相应 结果 。 由 表 1 及 图 5 和 图 6 知 ， 本 文 设计 的 非 线性 控制 涡 和 参数 白 适应 更 新 律 
对 有 界 噪声 具有 很 强 的 抗 于 扰 能 力 。 
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e(t) 
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a " 
25 30 
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a4(t), b,(0, c.t) 
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图 5: 系统 (9) 与 系统 (10) 的 改进 投影 同步 误差 ”图 6: 有 界 噪声 干扰 下 三 个 未 知 参 数 a, b, ctf eit 4l 


表 1: 系统 (9) 与 系统 (10) 最 大 同步 误差 及 对 未 知 参 数 估 计 的 最 大 相对 误差 


D absel 

1.2 0.0516 

2.0 0.0879 

4.0 0.1521 

6.0 0.2625 

8.0 0.3476 
5 fi 


absez 
0.0076 
0.0122 
0.0227 
0.0370 
0.0488 


abses 

0.1622 
0.3091 
0.5163 
0.9463 
1.1647 


E1 
0.7290% 
1.2526% 
2.3660% 
3.7263% 
4.9866% 


E2 
0.0587% 
0.0941% 
0.1737% 
0.2755% 
0.3968% 


E3 
0.0384% 
0.0333% 
0.0908% 
0.1132% 
0.0365% 


AX d 3e dad RE, E ih ER PET hA AC DE SUNL. WAT CORN 
混沌 系统 与 Genesio- Tesi 混沌 系统 的 改进 投影 同步 以 及 对 新 混沌 系统 未 知 参 数 估计 问题 。 
用 Lyapunov 稳 定性 理论 证 明了 Genesio-Tesi 混沌 系统 与 新 混沌 系统 在 改进 投影 同步 意义 下 误 
差 系 统 原 点 的 渐进 稳定 性 。 数 值 仿真 表明 本 文 设计 的 非 线 性 控制 器 和 参数 自 适应 更 新 律 能 够 快 
速 、 有 效 实现 新 混沌 系统 和 Genesio-Tesi 混沌 系统 改进 投影 同步 ， 准 确 估 计 出 新 混沌 系统 的 三 


个 未 向 | 参数 昌 对 有 界 品 声 干扰 具有 很 强 的 鲁 棒 性 。 
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Adaptive Modified Projective Synchronization and Parameter 
Estimation for Chaotic Systems with Uncertain Parameters 


LI Zhan-guo, XU Wei 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072) 


Abstract: A modified projective synchronization between a new chaotic system and a Genesio- Tesi 
chaotic system is investigated, and the unknown parameters of the new chaotic system are also esti- 
mated. The nonlinear controller and the parameter adaptive updating law are successfully designed 
according to the active control principle. Moreover, a strict proof of the asymptotic stability of the 
origin for the error system is derived by the Lyapunov stability theory. The effectiveness and the ro- 
bustness under a bounded noises disturbance of the proposed nonlinear controller and the parameters 
adaptive updating law are validated through numerical simulations. 

Keywords: new chaotic system; Genesio-Tesi chaotic system; modified projective synchronization; 


adaptive law; bounded noise 
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